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Abstract: 

There are two possible methods for the inclusion of the chemical potential in a relativistic bosonic 
^ ' field theory. The most popular method has recently been criticised by some authors, so much so as 

to require a re-analysis of the entire problem. I here present a new approach to the inclusion of 
the chemical potential that uses the second method and zeta-function regular is ation techniques. I 
first apply it to the non-interacting field, obtaining an expression for the effective potential which is 
formally coincident with the well-known one from the first method. My approach, however, seems 
f*"*. ■ here mathematically more rigorous since it excludes the presence of the multiplicative anomaly. I 

then obtain a new expression for the one loop effective potential for the interacting field. This is 
regularised, presents less cryptic excitation energies, shows Goldstone excitations, and in the non- 
fi |j interacting limit gives the expected expression. The high-temperature expansion is calculated easily 

^ ■ and with clarity, and confirms the results of the first method for the critical temperature. 



PACS numbers: ll.10.Wx, 05.30.Jp, 11.30.Qc 



1 Introduction 



The calculation of the one loop effective potential for a relativistic self-interacting bosonic field in the 
presence of a net global charge has been the subject of interest for a long time. The first calculations for 
the inclusion of the chemical potential in a relativistic field theory were done by Kapusta Q Q and by 
Haber and Weldon ||, in 1981. Kapusta calculated the partition function for a non-interacting charged 
scalar field, showing the presence of relativistic Bose-Einstein condensation. He also showed, using an 
approximation, that in the self-interacting case the chemical potential can modify the phase transition. 
Haber and Weldon confirmed these results and calculated, amongst other things, the high-temperature 
expansion of the thermodynamic potential Q |5| . 

These calculations, in the free field case, were later also developed using zeta-function and heat-kernel 
regularisation techniques || [jjj || and extended to curved spaces || Q ■ 



It was only recently though, with the work of Benson, Bernstein and Dodelson 12 |13|, that an exact 
expression was obtained for the one loop effective potential for a self-interacting charged scalar field. 
These authors included the chemical potential with the method, already used by Kapusta and by Haber 
and Weldon, and which I will call Method I, of inserting in the functional integrals that determine the 
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thermodynamic partition function, an effective Hamiltonian H — [iQ, where H is the microscopic Hamil- 
tonian, Q the charge operator and fj, the chemical potential. Having calculated the effective potential, 
they demonstrated, in the high-temperature limit, how the chemical potential noticeably changes the 
phase structure of the system, and confirmed the interpretation of spontaneous symmetry breaking as 
Bose-Einstein condensation. 

This remarkable work, however, has some unclear aspects. In particular, it is the approach itself 
for including the chemical potential that is criticised by some authors. Evans' analysis 1 14 shows how 
Method I is "unphysical" in the sense that it is not clear in the distinction between micro and macroscopic 
conditions, between evolution and boundary conditions and in the definition of the energies. There is 
also the problem of the definition of the charge operator in the case of broken symmetry. According 
to the considerations of Chaichian, Montonen, Perez Rojas and others, presented in a recent series of 
articles Jll| (l?]] , it is not legitimate to use, in the time-evolution operator, the charge of a symmetry 
that can be broken, because in these conditions the charge is not defined Jlq]|p^| . Benson, Bernstein 
and Dodelson's effective potential is then obtained by formal manipulation of divergent terms and is also 
not clear in defining the energies of the particles. For this latter reason Benson, Bernstein and Dodelson 
can only study the symmetry breaking in the high-temperature limit which is obtained in a decidedly 
complicated way. Method I, moreover, as I will show later, in the case of interacting field does not 
easily lend itself to regularisation with the zeta function and, in literature, it is necessary to resort to 
a heat-kernel expansion (20|. There remain also some doubts raised by Chaichian and colleagues jl7j 
concerning the renormalization in the presence of chemical potential carried out by Benson, Bernstein 
and Dodelson. These could be resolved by studying the Seeley-De Witt coefficients once an approach 
with the zeta function is possible. 

Besides the method of including the chemical potential used by Benson, Bernstein and Dodelson and 
those who preceded them, there exists another one which, as Evans recently observed |l4| , distinguishes in 
a much clearer way between microscopic and macroscopic quantities, keeping separate in the formalism 
the interaction Hamiltonian H and the chemical potential. Very little has been calculated with this 
method pl"|p2|, and it has never been used to calculate the effective potential. Following Evans' example 
I will refer to this alternative as Method II. 

From what was said regarding Method I it is clear that it is necessary to develop an alternative 
approach to the calculation of the effective potential, an approach which does not suffer from the same 
ambiguities of the existing one, that avoids, for example, the need to manipulate infinities, that permits 
the use of regularisation techniques, bypasses the problem of the charge operator and provides less cryptic 
excitation energies. The development of such an approach goes beyond being interesting, it is fundamental 
for verifying whether it is actually possible to include the chemical potential in an interacting field theory. 

The aim of this work is the critical comparison of the two methods, the analysis of the unclear aspects 
or of the difficulties in the practical use of both, but most of all the development of a new approach 
to the calculation of the effective potential, based on Method II and the zeta-function and heat-kernel 
regularisation techniques. 

In this work I will make extensive use of zeta-function regularisation, a technique introduced by 
Dowker and Critchley |p3| and Hawking Q , which is particularly useful when dealing with the ultraviolet 
divergences arising in the calculation of the one loop effective potential . This is related to the small 
disturbances operator A (i.e. the operator extracted by the quadratic term in the expanded action) by: 

CO -, 

V eff = + lndet(A) (1.1) 

eit (3V 2(3V y ' v ' 

where /3 is the inverse of the equilibrium temperature of the system. This determinant, as a product 
of eigenvalues, is formally divergent. The zeta-function regularisation gives it a meaning in the sense of 
analytical continuation. The zeta function £(s|^4) for a small disturbances operator A can be explicitly 
expressed starting from the Mellin transform of the trace of the heat operator related to A: 

1 r°° 

C(s\A) = — t^Tre- tA dt (1.2) 
r ( s ) Jo 
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This function can be analytically continued to s = 0, so that its derivative is well defined 
shown that the regularised one loop effective potential reads: 

S° C(0\A) C(0\A) 2 

where I 2 is a renormalization scale. 

In section 2 I will briefly introduce the two methods of including the chemical potential. In section 
3 I will present a new approach to the calculation of the effective potential applying it to the case of a 
non-interacting field and showing how it produces the same formal expression for the effective potential 
as the more classical method, but in a mathematically more rigorous way. In section 4 I will show the 
propagators obtained with the two approaches and will then re-propose a comparison of the two methods 
done by Evans. In section 5 I will briefly present the standard approach to the interacting case, analysing 
some of its weaknesses. In section 6 I will finally develop my approach to this case, obtaining a new 
expression for the one loop effective potential and analysing various aspects of it. 



2 The two methods 



I am interested in the calculation of the one loop effective potential for a theory described by the micro- 
scopic Hamiltonian density: 

A . A 



(2.1) 



where (f> 2 = <f>\ + (\>\ and the sum over the repeated indices i = 1, 2 is understood. This Hamiltonian is 
invariant under a global 0(2) symmetry. Noether's theorem then guarantees the existence of a conserved 
current J„ and the correspondent charge. The latter is formally defined as the space integral of the 
current Jo and is a constant of motion, and also a generator of the symmetry. In our case the charge 
density can be easily derived as: 



Q = l7 r 2 - 4> 2 tti 



(2.2) 



We therefore have a quantum operator Q associated with an additive quantity of the system and can 
define the grandcanonical partition function: 



Z = Tr{e-^ H ~^) = I d<j> «t>\ 



e 0(H-»Q)u> 



(2.3) 



It is now a matter of finding useful expressions for the calculation of this partition function, using 
functional integrals. 
In general: 



N 



[dcj>i][d(j>2\e 



(2.4) 



BC 



where BC are boundary conditions on the fields and S the classical action that is obtained from the 
chosen Hamiltonian. 
For Method I: 



S = 



dr / d A x 



1 



[(drfa) 2 + (drfo) 2 + (V^x) 2 + (V0 2 ) 2 + m 2 



A 

IT 



+ 



7 A 1 



(2.5) 



I have emphasised the term in /i 2 , which appears only due to the presence of the charge density Q in 
the effective Hamiltonian. This is a very important term for the development of the discussion. We can 
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get an idea of this by observing that this action presents an "effective mass" term m 2 — /i 2 . We expect 
therefore that this new contribution itself gives symmetry breaking when /i 2 is greater than m 2 , even for 
positive m 2 . 

The boundary conditions on the fields are the well-known periodic ones: 

0(t)=0(t + /3) (2.6) 

This approach does of course present some advantages. Including the chemical potential in the 
Hamiltonian, other than being intuitive, allows us to have an analogous theory to that without net 
charge and the periodic boundary conditions simplify the calculations. It is therefore understandable 
that this approach is more widely used, if not the only one used, for practical calculations Q [g || || Q . 

For the less well-known Method II instead, one obtains: 





dr I df'x 



{drfaf + {d T <hf + (V0l)" + (V4>2? + m 2 ct> 2 



(2.7) 



Now the chemical potential is not present in the action and only the true microscopic Lagrangian appears. 

Instead of the usual periodic conditions, we have pseudoperiodic conditions on the fields, that contain 
the chemical potential [^l| |2^] : 

0(r) = e^0(r + 0) (2.8) 

These boundary conditions are less convenient than the periodic ones and it is probably this that has 
prevented a greater use of this method but, as I will show later, this is balanced by much simpler small 
disturbances operators. 



3 The non-interacting field 

The case of a non-interacting field has been widely studied with Method I. It remains interesting, though, 
for understanding where the characteristic problems of the interacting case arise. It is also an important 
test bench for any attempt at a new approach to the problem. 

With regards to Method I, it can be easily seen from 2.5 that with the non-interacting action is 
associated an operator Aj\ 



a 2 - V 2 + m 2 - [i 2 -2ifid T 



2ind T -d 2 - V 2 + to 2 - /i 2 



(3.1) 



This is a matrix valued elliptic differential operator, but not self-adjoint, due to the presence of the 
chemical potential. 

I am interested in the effective potential, which is linked to this operator by: 

Veff= 2jW lndet(Al) (3 - 2) 
As Ai is matrix valued we formally have: 

Veff = ^lndet[(-a 2 -V 2 + m 2 - M 2 ) 2 -4 M 2 a 2 ] = 

' lndet{[(id T - ifi) 2 - V 2 + m 2 ][(id T + i/x) 2 - V 2 + m 2 }} = 



2f3V 



= ^^lndet[(ia T -z^) 2 -V 2 +™ 2 ] + -J-lndct[(ic) T +^) 2 -V 2 +m 2 ] (3.3) 

zpv zpv 

It is now easy to find the eigenvalues of these operators as: 

(ui n ± ifi) 2 + k 2 + m 2 (3.4) 



4 



with oj n = ^S-n. 

With these eigenvalues: 



Tre 



-tAi 



V 



d 3 k 

(2^)3 



-t[(u„+i/i) 2 +fc 2 +m 



-t[(u 



-i/j,y+k 2 +m' ! 



(3.5) 



From this trace of the heat kernel the zeta function can be easily calculated and hence also the effective 
potential. 

Infact attention has only recently been drawn p5[ to the fact that two zeta-function regularised 
determinants det(B) e det(C) do not satisfy, in general, a relationship such as lndet(_BC) = lndet(-B) + 
lndet(C) (3.3) and to assume its validity (as has always been done up to now) can cause a quantity called 
multiplicative anomaly to be neglected. This quantity is given by: 



\nF(BC) = hidet(BC) - lndet(B) - lndet(C) 



(3.6) 



Such a quantity may be found to be physically relevant, but it seems not easily computable, especially 
for not self-adjoint operators like those with which we are working. 

We will now look at how it is possible to calculate the non-interacting effective potential with a new 
approach that uses Met hod II and the zeta-function regularisation techniques. 

To the action of 2.7, for A = 0, is associated an operator Ajj: 



-dl 



v 2 
o 





-dl - V 2 



(3.7) 



This is a matrix valued elliptic differential operator, not self-adjoint. An important characteristic for us, 
with respect to the corresponding operator calculated with Method I, is that here there are no differential 
terms on the off-diagonal. As we will see later, this will greatly simplify the calculations with the zeta 
function in the self-interacting case. 

My approach also has another remarkable advantage with respect to Method I: the operator that has 
just been calculated is made up, on its diagonal, of two identical opcratorial terms and therefore does 
not present any multiplicative anomaly. 

Even if I have an operator which is not self-adjoint, but is elliptic, I can still define the logdet(A//) 
via the zeta function and I can calculate it explicitly as I know its spectrum. Even though the chemical 
potential does not appear explicitly in the operator, because of the boundary conditions it appears in the 
couples of eigenvalues: 



(u> n + ifl) 2 + k 2 + TO 2 

The trace of the heat kernel will therefore be 

d 3 k 



(3.8) 



Tre 



-tAi 



V 



2V 



(2tt) 3 
d 3 k 



-t[(w„-MM) 2 + k 2 +m 2 ] , -t[(w rl +i/i) 2 +fe 2 +m 2 ] 



(27T) 



e -t[(w„+iM) +k-+m 2 ] 



3 



(3.9) 



It is therefore possible to continue with the standard methods of zeta-function regularisation and obtain 
the well-known: 



1 



eff 



16tt 2 2 
I f d 3 k 

1, 



(27T) 



- ln(m 2 £ 2 ) + - 



ln(1 _ e-W-rt) + 1 

P 



d 3 k 
(2^P 



ln(l 



= _ /3 ( v /fe2 +m 2 +/ _ I ) 



(3.10) 
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This expression coincides with the analogous ones calculated with Method I but, unlike those, it does not 
require any eventual additional terms due to the anomaly. 

It is now easy to see how the two approaches anyhow coincide except for the possible multiplicative 
anomaly. 

The eigenvalues of the two operators are different, but what is important for us is the trace of the 
heat kernel. 

With Method I: 

d 3 k 



Tre~ tAl = V 



j ( t)3 £ [ e " 



but: 



(uj- n - ifi) 2 + k 2 + 



2tt 

(~—r n — *M) 2 + k 2 + m — (uj n + i/i) 2 + k 2 + 

P 



m 



and the summation is for all the integer numbers. 

Hence, also for Method I, the trace of the heat kernel (except for the anomaly) is: 



Tr e 



-tAr 



2V [ rf3fc y^ c -t[(^+»M) 2 +fc 2 +r» 2 1 = Trt 
J (2tt) 3 ^ 



-tAi 



(3.11) 



(3.12) 



(3.13) 



Once the trace of the heat kernel has been calculated, via a Mellin transform one obtains the zeta 
function and from this the effective potential. Because this is the only object from which we can obtain 
physical information about the system, we can see how the two approaches, in the charged free case, give 
coinciding results (except for the multiplicative anomaly) , my approach being more rigorous and reliable 
if the anomaly is found to be influential in this case. 



4 Propagators 

An analysis of what I have introduced so far can show how Method I is lacking at least in physical 
clarity. As we have seen, infact, Method II keeps that which determines the microscopic interactions, the 
Hamiltonian of the zero-temperature theory, well separate from that which determines the "macroscopic" 
boundary conditions, the chemical potential. We have, infact, taken as the complex time-evolution 
operator the term e~ rH and only with this have we established the relationship between partition function 
and functional integrals. The chemical potential, instead, only appears in the boundary conditions. 
Method I, on the contrary, inserts the macroscopic information together with the microscopic ones directly 
into an effective Hamiltonian which it then uses both for the thermodynamics as well as for the microscopic 
evolution of the system, with a complex time-evolution operator e~ T ( H ~ fJ - < 3) . 

The problem is not just about "aesthetic" differences or, as we also shall see, of greater or lesser ease 
of calculations. In the near-total absence of material in literature, a recent comparison was made by 
Evans Q. Evans gives substance to these observations with a comparison of the propagators, calculated 
with both methods. 

It is possible to calculate the thermal propagator (euclidean signature) also with my approach by 
integrating on the parameter t in the heat kernel associated with the operator An. 

I here show only the result, which is analogous with that already present in literature and calculated 
with more classical methods |Q |2l| . 

f dk f d 3 k e £3+lkaT 
G{x,0) — / — — 



(4.1) 



2tt J (2tt) 3 k 2 +uj2 



d 3 k „ j& f 1 e"— 1 



re 



(2tt) 3 ^ e /3(ofe-M) - 1 2uo R e^^+c) - 1 2ufe 



G 



where ujfi = k 2 + m 2 . 

The first term of the propagator is the zero-temperature contribution which shows us the poles. 
The other two are clearly interpretable as the thermal contributions. There appear, infact, Bose-Einstein 
distributions for particles and anti-particles. Most of all, it is these terms that contain all the dependencies 
from the macroscopic parameters which, as one would expect, do not appear in the poles. We have infact 
two poles, in the same position in which we would find them if in the absence of chemical potential: 

h = ±iufe (4.2) 

The poles, therefore, do not depend on the macroscopic physics but only on the form of the microscopic 
Hamiltonian, the true physical Hamiltonian. 

Consider now the propagator for Method I. It can be easily seen that it has two poles in: 

fc = i(±ofc - /*) (4.3) 

This is quite a peculiar result as the poles depend on the chemical potential; we therefore have the 
presence of a parameter characteristic of "many-body" in a microscopic term. In his work Evans shows 
how this paradox can be recomposed by studying the relationship between the fields, in the Heisenberg 
representation, used in the two methods. These are not the same due to the different definitions of the 
evolution operator. The relationship between two Green functions is such that the poles are related by a 
shift of \x in the energy. 

It is important to remember at this time that the two methods are defined by starting from the same 
partition function and therefore should provide the same physical information about the system, even if 
they can differ in the way they do this. Method I, though, requires a much greater attention. The energies 
of the various particles are not measured with respect to a common zero but are traslated, through the 
choice of H — \iQ, by an amount equal to q/j, with respect to the standard zero of relativistic field theory 
and of Method II itself. This could lead to some confusion, also when, for example, we are dealing with 
many charges and symmetry groups, as in the standard model. 



5 Interacting field: standard approach 

In the background field tp approximation, the small disturbances operator for Method I is: 

/ -d 2 - V 2 + to 2 - [i 2 + \ V \ + §<p| -2i/z9 T + !<p l¥ > 2 

V 2i[id T + §<pi<p 2 -d 2 -V 2 +m 2 ~(i 2 + ±tp 2 + ±tp 2 

It is a matrix valued, non-self-adjoint, elliptic differential operator. 

If we compare this with the analogous expression for the free field, we notice immediately the presence 
of the terms in A on the off-diagonal. It is actually the presence of \i and A together in that position that 
complicates things. The eigenvalues of the operators associated with lndet(v4j) are infact: 

col + k 2 + m 2 - M 2 + ^ 2 ± yj - 4 M 2 W 2 (5.2) 

If we did not have the interaction the root would disappear and the same would occur if there was no 
net charge. We know, infact, that with Method I problems arise only for the interacting field. In this 
state, though, the eigenvalues cannot be manipulated easily. The standard way of proceeding is to write 
the trace of the heat kernel and use the Poisson summation formula. Here, though, the root is a great 
obstacle and one does not obtain much. Not even an approximation for high-temperatures made directly 
on the eigenvalues produces useful results. In literature (even if they are calculations on curved spaces) 
usually at this point one has to limit oneself to the expansion of the heat kernel and then proceed by 
approximation p(| . 
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The one loop effective potential calculated by Benson, Bernstein and Dodelson through formal ma- 
nipulation of infinities is anyway: 

1 , 2 2 A 4 , f d 3 k E + (k)+E_(k) 
Veff = ^(m -n)<p +jj<p 

1 f d 

+ tJ-2 

with: 



2 V /r ' 4! r ' J (2tt) 3 2 

1 ' ^Ml-e-^ (fc) ) + i / |* Ml-e-^) (5.3) 



£ ± (fc) 2 = fc 2 + m 2 + M 2 + ^ 2 ± J + 4 M 2(fc2 + m 2 + * 2) (5.4) 

3 V 36 3 

Benson, Bernstein and Dodelson themselves study only the various limits of this result. The non- 
interacting limit coincides with that expected and, in the limit fi tending to zero with broken symmetry, 
one finds Goldstone excitations. The high-temperature limit turns out to be much more complicated. 



6 Interacting field: a new approach 



The use of Method II for the interacting field is not as easy as for the free field. Due to the new 
pseudoperiodic conditions on the fields infact, it is not possible to proceed with the usual expansion with 
respect to a constant <f>. The background field (f> must be dependent on the imaginary time. It is infact 
possible to study the relationship between the fields in Method I and in Method II |TJ] and, keeping in 
mind the different definitions of the evolution operator in the two cases, one obtains: 



4>ii{t) 



Mr) 



(6.1) 



where I have labeled the proper fields of each method. In particular, when one studies a system with 
symmetry breaking which has a background field <fi^ 0: 



hi W 



h (r) 



(6.2) 



with ip constant as the background field in Method I is constant. From now on I will omit the index II. 

Let us see what this dependency implies in the context of an expansion of the action. If I expand the 
action up to the second order in the fields with respect to <f>, I obtain: 




Vi(x) + | / d 4 xd i y r]i(x) 



S 2 S 



5<j)i{x)84>j(y) 



With the necessary functional derivatives and remembering that —V 2 ^ (r) = we obtain that: 



and: 



5 X = 



5 2 = ; 



dr / era; 



\ 4> t (r)(-9 T 2 ) & (r) + \m> f (r) + A £ ( T ) 



<It j <f.r j -d 2 T fa (r) +m 2 fa (r) + ~(<g (r)+ <g (r)) fa (r) j ,/,(,•) 



dr Id x rj t (x) A it j r]j(x) 



(6.3) 

(6.4) 
(6.5) 

(6.6) 
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where Aij is the matrix that represents the small disturbances operator A: 



-i) 2 -Y 2 + m 2 + ^fa+^fa 

A A A _Q2 Y72 i ™2 | A ~" _l .. :'l 



2 _ Y72 j_ ™2 j_ A A _L A A A... 

2 ,-2 1 (6-7) 



\ 3 VI V2 ~ u t — v "6" VI ^ 2 ^2 / 

Let us observe first the operator A. It is a matrix valued elliptic differential operator. It is not self-adjoint, 
not because of its form, but due to the boundary conditions, whereas the analogous Aj for Method I 
had "good" boundary conditions but the presence of the chemical potential did not render it self-adjoint 
anyway. This operator has a great merit compared to that of Method I: the differential terms which 
caused so many problems do not appear on the off-diagonal any more. It would seem therefore that this 
approach greatly simplifies the calculations and the eventual results. Unfortunately there is a price to 
pay for this: the term (j) (r) is also in the small disturbances operator. 

With regards to the term S°, we would expect to obtain from it the constant classical potential, 
something which does not happen if the <fi depends on time r. 

Also the term S 1 causes us some problems. I here show fully the equations of motion: 

-d 2 T + m 2 + ^{fa (r)+ fa (r))) fa (r) = 

< + m 2 + ^{fa (r)+ fa (r))) fa (r) = (6.8) 

These are the conditions necessary for the term S* 1 to be cancelled and so that we can proceed as usual. 
It can be easily seen, though, that our cfr (r) = e^^cp, if it is not constant (i.e. jj, = 0) or null, cannot be 
a solution to these equations. We would infact have: 

V + ™ 2 + \{fa (t)+ fa (r))) e-^ifi = (6.9) 

and therefore, unless is zero: 

(fa (t)+ fa (r)) = y(m 2 - m 2 ) = constant (6.10) 
A 

The only exact solution to the equations of motion (for [i ^ 0), therefore, is obtained for a null background 
field, i.e. when the simmetry is unbroken. 

I was able to overcome the obstacles described above by resorting to a high-temperature expansion, 
which is anyway what we are most interested in as it is in that range that we want to study phase 
transition phenomena. As I showed before: 

4>= e-" T ip = (1 - iit + ^ 2 t 2 + . . .)ip (6.11) 

For high temperatures (3 is small and therefore r is also. We could therefore limit ourselves to substituting 
4> with its zero order approximation in r, ip. This would seem, at first sight, a good solution to our 
problems with the equations of motion. Doing this, however, we eliminate very important physical 
information. The tree term of the effective potential would not be dependant on the chemical potential 

-2 

anymore. Whereas with Method I the chemical potential is explicitly present in S° as ^/j, 2 <j> , here it is 

introduced due to the dynamic term —d 2 4> and therefore would not appear if we fix <f>= (p. 

Not even stopping at the second order in the expansion of the exponential is of any use. If we obtain 

the desired term in the dynamic part, there remains, however, the time dependence in the terms of the 
, -2 

type m <p . 

Let us take instead <j> (r) and substitute it in the equations of motion. We obtain: 

(- M Vi + O(r)) + m 2 (l — jir + i^r 2 + • • + ^4? (r) fa (r) = (6.12) 
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I limit myself to taking the zero order in t of the equations of motion: 
A 



H 2 ifi + m 2 ip l + -tp 2 <pi 
6 







(6.13) 



Now a non-zero tp can be a solution. I have not, therefore, done an expansion for high temperatures of 
the background field, but of the equations of motion themselves, avoiding in this way the loss of relevant 
physical terms and showing how in the high-temperature approximation the term S 1 is negligible. The 
same can be done with the zero order term of the action: 







dr I d x 



0V[- 



1 2 2 i 1 2 2 

"2^ f + 2 m V 



1 2 2 , ^ 2 2 i ^ 4l 

■-tup +-m <p + 



A 4 
4!^ 



and with the small disturbances operator: 



-d 2 - V 2 



2^1 



1^2 



-ffi - v 2 



(6.14) 



(6.15) 



I will proceed therefore in the calculation in the one loop approximation of the effective potential 
using this operator, which, unlike Method I, allows me to exploit the zeta-function techniques. For the 
form of the operator we cannot here exclude a priori the presence of multiplicative anomaly. 

We must remember that, in the case of broken symmetry, it is only a high-temperature approximation. 
It is, instead, an exact result when the symmetry is not broken. 

As expected, the eigenvalues of the operators associated with A, calculated in an analogous manner 
to 3.3, have a particularly simple form: 



K + + k 2 +m 2 + ^ip 2 ± ^tp 2 
3 6 



(6.16) 



I will avoid here, for brevity, the explicit calculations. These can be found in the appendix. The one 
loop effective potential, calculated using the approach for the inclusion of the chemical potential that I 
developed, results finally as: 
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(6.17) 



where M 2 = m 2 

It was obtained using zeta-function techniques and, as such, is regularised. It is to be noted that the 
excitation energies are all particularly clear. 

This result naturally includes within itself also the case of non-interacting field. In the limit of A 
tending to zero we obtain infact, for the effective potential, the exact expression: 
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(6.18) 



This expression coincides with that calculable with my approach starting directly from a free field, in 
other words using the 3.9 for the |3.10[ As such it coincides, except for the possible multiplicative anomaly, 
to the expression calculated in the classical manner and widely studied in literature. 

We know that this is also an exact calculation for the effective potential when the symmetry is 
unbroken. If (p— ip = 0, then M 2 = m 2 and we easily obtain a coincident expression to that of Benson, 
Bernstein and Dodelson's analogous limit. 

If we are in conditions of broken symmetry, this expression cannot be seen as anything other than a 
high-temperature approximation. 

Let us consider the case, with m 2 < 0, of symmetry spontaneously broken even at zero temperature. 
Now, the value of tp that minimises the tree- level potential is necessarily such that: 
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(6.19) 



Substituing it in the thermal excitation energies of 6.17, we get: 

Ei(k) -- 



X 



ip 2 ± fi = \Jk 2 — 2m 2 + 3/x 2 ± [I 



E 2 (k) 



k 2 



ip 2 ± fx = \/k 2 + [i 2 ± [i 



Now taking the limit for /i tending to zero: 
E ± (k) = ^k 2 + 2\m 2 \ E 2 {k) = \ 



(6.20) 



(6.21) 



As we can see from 3.2, for /i — the background field is constant and hence we return to an exact result. 
We can therefore interpret these expressions as the energies corresponding to the expected Goldstone 
boson and a particle of mass y / 2\m 2 \. 

I would like to point out that the expression obtained by Benson, Bernstein and Dodelson, though 
it may have been exact, in its complexity does not allow us to extract much more information about 
the system than those obtained with my result. The authors themselves, must limit themselves to a 
high-temperature expansion when studying the phase transition, the influence of the chemical potential 
and when giving an evaluation of the critical temperature. Moreover, their method of calculating this 
approximation is rather complicated and requires, for example, the analysis of two different regimes 
of charge density. With the expression I have just derived, instead, obtaining the high-temperature 
expansion is nearly immediate. 

Let us observe, infact, that the thermal terms of 3.17 have the same form as those of the free case: 
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The high-temperature expansion of this expression was first calculated by Haber and Weldon || 
well-known and now easily obtainable using a Mellin-Barnes representation of the zeta function 
result is, in the leading terms in T: 
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If I use the evident analogy between my result and the free case I obtain: 
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This result coincides, disregarding 0(T) orders in the thermal part, with that obtained by Benson, 
Bernstein and Dodelson. The terms in T 4 , T 2 arc dominant and one usually limits oneself to these for 
studying the phase transition. The difference between Benson, Bernstein and Dodelson's result and mine 
appears only from the term in T of the thermal contribution onwards, but this is reasonable as they are 
both approximate calculations. 

Benson, Bernstein and Dodelson's estimations for the transition temperature and the interpretation 
of the spontaneous symmetry breaking as Bose-Einstein condensation are therefore confirmed by my 
approach here. 



7 Conclusions 

I have here analysed the two different methods for the inclusion of the chemical potential. I showed 
the existence of some weaknesses of the standard approach, demonstrating how the development of a 
new alternative approach to the problem is not only interesting but necessary for verifying the actual 
possibility itself of including the chemical potential in an interacting theory. 

Such an approach was developed using Method II and the zeta-function and heat-kernel regularisation 
techniques . 

Applied to the case of a non-interacting field I showed that it was equivalent to the well-known and 
widely used method if the multiplicative anomaly, which could complicate the result of Method I, was 
disregarded. It, instead, does not appear in my approach which seems here mathematically more rigorous. 

I then obtained the propagator and compared it with that calculated with Method I, showing how 
the new approach distinguishes more clearly between microscopic dynamics and boundary conditions. 

I then applied the new approach to the more complicated interacting case. A comparison between 
the small disturbances operator in the two approaches showed that the new method was more suitable 
for a regularised calculation. I obtained then a new expression for the one loop effective potential for the 
interacting field. It is the first expression for this effective potential obtained using Method II, a method 
that, as I highlighted and according to other authors, does not soffer from some of the ambiguities of 
Method I. The new expression is regularised but, in the case of broken symmetry, is a high-temperature 
approximation. The corresponding excitation energies are also less cryptic than those obtained by Benson, 
Bernstein and Dodelson with Method I. 

I then studied some relevant limits. In the limit of non-interacting field I again obtained the expected 
expression, and likewise in the case of unbroken symmetry. In the case of symmetry breaking with negative 
mass squared, the expected Goldstone excitations appear. Above all, the high-temperature expansion of 
my expression turns out to have been obtained in a much simpler and clearer way than that of Benson 
Bernstein and Dodelson. My expression coincides with it, however, in the leading terms in T, confirming, 
therefore, the results obtained by these authors for the symmetry breaking at high temperatures. 

This new approach could be further developed. For example, an analysis of the Seeley-De Witt 
coefficients could resolve certain unclear aspects of the regularisation in the presence of chemical potential, 
it could be applied to curved spaces and the simpler excitation energies could be of great help in problems 
of dynamics of out of equilibrium systems, where knowing the "masses" of the particles in play is crucial 
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Appendix: Explicit calculations 

To make the discussion more compact I here present the necessary calculations for obtaining the inter- 
acting one loop effective potential using my approach. Although it could appear complicated, these are 
standard calculations of the zeta function technique. 
I define: 

Mf = m 2 + ^ 2 Mf = m 2 + ^ 2 (A.l) 

Substituting in |6. 1€ , the trace of the heat kernel is therefore: 

T re - tA = V [ ^[c- t[( "" +tM)2+fc2+Ml21 I c -*f(^+'M) 2 +fc 2 +Mjlj ( A 2 ) 

J ( 27r ) 3 „ 

Using Poisson's summation formula [^| : 

E f( k + l)= E / f(x)e 2niHx - q) d N x (A.3) 
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The zeta function will then be: 
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(A.5) 



I have distinguished the vacuum term (n = 0) from the thermal contributions and have defined uM = 
k 2 + M 2 . As the two terms relative to the two eigenvalues differ only in "mass" , I will only explicitly 
indicate one of them. 
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Let us consider first the vacuum term. Remembering the definition of the function T(s): 
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In the last step I used T(3/2) = y/n/2 and the identity: 
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This vacuum contribution is analytic in zero. I can therefore easily differentiate and obtain: 
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I take back two of the terms of the thermal contribution: 
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I use the representation |2j| of McDonald's functions: 
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With z = Pnufe e v — 1/2 — s: 

It is an analytic function. I obtain, therefore, for one of the two cases: 
Remembering another useful representation of McDonald's functions [E 
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we then pass easily to: 
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And from here to: 
, d 3 k 

V 



ln(l 



The condition for convergence of the sum is here turned into: 
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The term with the negative chemical potential, proceeding analogously, results in: 
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(A.17) 



with the condition: —fi < M\. 

We must remember, though, that we have two terms in M2 that give analogous contributions. By 
definition, M 2 < Mi and therefore the condition for summation this time requires that the chemical 
potential is: 



< \Jm 2 + ^cp 2 
I therefore found that: 
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To determine the effective potential we need also £(0|A) ln(^ 2 ). As all the elements of the thermal 
contribution are of the form f(s)/T(s), these do not contribute as f(s) is an alyti c and therefore finite in 
zero, whereas T(s) has a pole in zero. Therefore, remembering that £o(s|^4) ( A. 6 ) is analytic in zero: 



C(0|^)ln(£ 2 ) = Co(0|A)ln(£ 2 



(3V Mf 
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(A.20) 



Combining these contributions we obtain my expression 6.17. 
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